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Abstract 

The evolution with energy scale of the partonic logarithmic scaled energy spectra is investigated in the 
framework of the modified leading logarithmic approximation (MLLA). The behaviour of the higher order 
moments is compared to a number of analytic predictions and e + e~ data. 
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1 Introduction 



Quantum Chromodynamics (QCD) currently does not allow direct calculations of the hadronic final state 
observables. To make predictions of the final state it is necessary to model the transition from partons to 
hadrons. One approach that has been successful in describing the general features of the inclusive energy 
spectra, in both e + e~ annihilation and deep inelastic scattering experiments is the modified leading log 
approximation (MLLA) using local parton hadron duality (LPHD) as the method for relating the MLLA 
partonic predictions to the hadronic observables. 

The production of hadrons in hard scattering processes is controlled by the underlying partonic behaviour. 
At small values of momentum fraction (of the outgoing parton from the original hard scatter) this parton 
behaviour, often referred to as a parton shower, is dominated by gluon bremsstrahlung. The branching processes 
q — > qg and g — > gg (double logarithmic processes) in addition to g — > qq (single logarithmic process) give rise 
to this parton shower. MLLA accounts for both these double and single logarithmic effects in the evolution 
equations The perturbative properties of partonic distributions have been calculated in the framework of 
the MLLA. They are governed by two free parameters: a running strong coupling, defined by a QCD scale A, 
and an energy cut-off, Qq, below which the parton evolution is truncated. 

Using LPHD, the non-perturbative effects of particle distributions are reduced to a simple factor of normal- 
isation that relates the hadronic distributions to the partonic ones. At large enough energies, away from the 
influence of Qq, this 'hadronisation' factor should be independent of the energy scale [|| at which the spectra 
are being calculated. 

Various approaches, detailed below, have been taken to calculate the single particle spectra and their mo- 
ments within the MLLA framework. In this paper these theoretical approaches are compared as a function 
of energy scale, Q. Substantial difference are found in the predictions of the higher moments. The theoretical 
results, in particular the higher order moments, are also compared with e + e _ data over a range of centre of 
mass energies, Q = Ecm/2- It is found that the general characteristics are well described but that the theory 
is not in accord with all features of the data if it is assummed that A, Qo and the LPHD normalisation are 
energy independent. The studies reported here avoid the need to extrapolate into regions not experimentally 
measured. 



2 Single Particle Spectra 

Given a high energy parton which fragments via secondary partons into a jet of hadrons, the MLLA evolution 
equation allows the secondary parton spectra for the logarithmic scaled energy, £, to be calculated [Q. The 
variable £ is defined as \n(Q/E) = ln(l/x), where Q is the energy of the original parton and E is the energy of 
the secondary parton. The cut-off, Qq, limits the parton energy to E > k^ > Qq, where kx is the transverse 
energy of the decay products in the jet evolution. In order to reconstruct the £ distributions an inverse Mcllin 
transformation is performed 



where the integral runs parallel to the imaginary axis on the right of all singularities in the complex oj— plane, 



The Mellin-transformed distributions, D(u, Y, A), can be expressed Q in terms of confluent hypergeometric 
functions, <&, as 
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and A and B are defined as 



A = 4N c /buj, 



B = a/b, 
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where N c is the number of colours, a = ll-/V c /3+2n^/3iV ( ? , nt is the number of flavours and b = HN c /3+2nf /3. 
Throughout this paper the number of flavours is assumed to be 3. Equation ([l]) is then calculated using a 
numerical integration in the complex cj— plane. 

It is convenient to investigate the MLLA spectra in terms of moments. The cumulant moments of the £ 
distribution can be written as: 



K m = - 
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The analytic form of these cumulant moments have been calculated in B for the first four moments. This allows 
the normalised moment, (£ m ), to be calculated and hence the dispersion (a), skewness (s) and the kurtosis (k) 
to be constructed, see for example B 
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The so-called limiting spectrum is the case when Qq = A i.e. A = 0. For this case, Fong and Webber || 
have also calculated the behaviour of the moments of the £ spectra with energy scale, Q. They point out that 
the spectra can be represented close to the maximum of the distribution by a distorted Gaussian of the form: 
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where S = (£ — I) /(J and I is the mean. This expression, when fitted to the spectrum, allows the moments to be 
determined when the full spectrum is either unmeasured or uncalculable. 

The changes to the parton spectra for quark, q, and gluon, g, jets are Next-to-MLLA effects which result in 
additional terms in the integral shown in Eq. ([TJ) . These corrected spectra (at Qo = A) for the quark and gluon 
jet can be related to the limiting spectra M as follows 
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with 
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A,j. g . The mean is the same as 
7 a . These shifts are small with 



The effect on the moments is that the a, skewness and kurtosis for the quark/gluon distribution would be 
approximately the same as the limiting distribution at an effective energy Y + 
that of the limiting spectra at this effective energy but shifted in value by —A 
the limiting spectra being between the distribution for the quark and gluon jets. It should be noted that the 
predictions of Fong and Webber @ for the relative shifts between quark and gluon jets are slightly different 
from those discussed above. 

In the MLLA approach, the partons are assumed massless so the scaled energy and momentum spectra are 
identical. Experimentally the scaled momentum distribution is usually measured and as the observed hadrons 
are massive the equivalence of the two spectra no longer holds. In || the assumption is made that the cut-off 
Qq can be related to the masses of hadrons. This allows the logarithmic scaled momentum distribution, £ p , to 
be written as 



1 dn h p h - 
~N~dip ~Eh 



(12) 



where 

^ = ^ y / Q 2 e~ 2 ^ + Ql ' 

the energy of a hadron with a momentum ph is .E/j = \Jp\ + Q 2 - Limiting momentum spectra based on massless 
partons and massive partons will be referred to as MLLA-0 and MLLA-M spectra, respectively. 



3 Behaviour of Theoretical Spectra 

In this section the general characteristics of the spectra are discussed in order to illustrate their behaviour as a 
function of the three variables Q, A, and Qq. 

Figure [|(a) shows the limiting energy spectra for three energies Q = 183. /2 (LEP II), 91.2/2 (LEP) and 
14.8/2 GeV (TASSO) with a A value of 250 MeV. As the energy scale increases, the parton multiplicity (the 
area under the curve) grows and the peak position shifts to the right i.e. the parton scaled energy spectra is 
softer. The cut-off on the right hand side of the plot corresponds to the truncation in the parton energy spectra 
at Q = A. 

Figure |l|(b) shows the scaled energy spectra at fixed energy scale Q = 91.2 GeV and A = 250 MeV but at 
two different values for the cut-off of the parton evolution, Qq = A and 2A. (The curves are cropped at £ = 1.0 
as below this value the numerical integration of the truncated spectra is unstable.) Truncating the cascade at 
higher values of Qq leads to a lower parton multiplicity with a harder energy spectra. 

Figure [l](c), again, shows the limiting energy spectra but at a fixed energy scale, Q — 91.2/2 GeV, with 
values of A = 50, 250 and 400 MeV. As A decreases, more partons are produced and their energy spectra is 
softer. In the case of the limiting spectra, this behaviour is dominated by the fact that an increase in A is 
accompanied by an increase in the cut-off Qq. 

Figure [l](d) shows the limiting scaled momentum spectra for the MLLA-0 (Eh = Ph) and MLLA-M case 
(Eh ^ Ph) at an Ecm = 91.2 GeV. The MLLA-0 spectra displays the usual truncation at large £ associated with 
the cut-off Qo, but the MLLA-M case does not. This is due to the fact the calculation is now regulated by the 
cut-off entering as a mass term in the expression of £ p (see Eq. 
introduction of a mass term has no effect on the MLLA calculation. As £ p increases (the momentum becomes 
smaller) the mass term begins to play an increasingly important role. From these arguments it can be seen that 
as Q decreases mass term has a more significant influence over a larger (fractional) range of the £ p spectra. 

Figure || shows the evolution of the mean, a, skewness and kurtosis of the limiting energy spectra with 
A = 250 MeV as a function of Q. Three different approaches to calculating the moments have been investigated: 

(i) the analytic expression of Fong and Webber for gluon jets (dashed line) and quark jets(dash-dotted line); 

(ii) the analytic calculation of Dokshitzer et al. Q| (dotted line); (iii) fitting a distorted Gaussian (Eq. [l(]) over 
±1(7 around the mean value (full line). The range of the fit is motivated by the phenomenological Q scope of 
validity. All approaches exhibit the same trends, in that all the cumulants increase as Q increase. The values of 
skewness and kurtosis are negative, tending towards zero as Q increases, i.e. the spectra are becoming more like 
a pure Gaussian. The mean and the a exhibit a very similar dependence on Q for all three approaches, though 
with different offsets. The skewness and kurtosis exhibit a different dependence on Q, but as Q approaches the 
asymptotic limit the predictions are converging. 

The differences between quark and gluon jets for the predictions of Fong and Webber, in Fig. |[ are small. 
A feature of the data is the marked difference between the analytic calculations of Dokshitzer et al. and the 
fit to the limiting spectra. Calculating the moments over the full range of the spectra gives the same result as 
the analytic calculations, not too surprisingly since both calculations use the same assumptions. Investigating a 
limited range of the spectra, as is done during the fit, produces marked differences, highlighting the sensitivities 
of these higher-order moments to the behaviour of the tails as well as the incorrect form of the distorted Gaussian 
away from the peak position. The difference between the analytic calculation of Fong and Webber and that 
of Dokshitzer et al. could well be due to the differences in the predictions at large x p , where the theoretical 
assumption made in both calculations are no longer valid. 

Figures [j] and |^ imply that A influences the position of the maximum and the width of the distribution of 
£. The effect of the relation between Qq and A is more complex; setting Qo = 2A lowers the position of the 
maximum for a given A. The maximum can be returned to the original position by reducing A but with a 
consequential increase of the width of the distribution relative to that found for the limiting spectrum. The 
effect of introducing mass terms into the momentum spectra is to increase the tail of the spectra at high £, thus 
increasing skewness. As Q increases, the influence of the mass term becomes less important. 



12 .) As the momentum — > Q, i.e. £ p — > 0.0, the 



4 Comparison with data 



The differential £ p distributions for e + e _ data and the models are shown in Figure g|. The limiting spectra 
were calculated with A = 250 MeV; this value was chosen to give the correct peak position of the spectra for 
an Ecm{= 2Q) at the mass of the Z°, Mz- For each energy, the theoretical spectra are normalised to the same 
maximum height as the data. 

The theoretical spectra follow the general trend of the data for both for both MLLA-0 and MLLA-M spectra; 
in detail there are discrepancies. At low energies, for the MLLA-M case, the skewness is larger (i.e. more positive 
) than that observed for the data; at high energies, the limiting spectra have a smaller width than the data. Also 
shown are the MLLA energy spectra for Qo = 2A. Here A = 50 MeV is required to get the correct position of 
the maximum for Ecm — Mz- The position of the maximum has approximately the correct energy dependence; 
however, the width of the theoretical spectrum is consistently greater than that of the data. 

To study the energy evolution of the MLLA theoretical spectra and the e + e~ data, fits around the peak 
position, using the distorted Gaussian (Eq. |l0|) have been made to determine the moments of the £ p distributions. 
In view of the statistical limitations, the range of the fit for the data was about three units around the peak 
(see Figure ||). The distorted Gaussian describes the data well over the fitted region with a \ 2 /dof of 1.0 or 
better. To permit direct comparison with the data, a similar fit range was used for the theoretical spectra. 
The distorted Gaussian gives a good description of the MLLA-M spectra. The description is less good for the 
MLLA-0 spectra but for all energies the Gaussian represents the model to better than 1%. 

The moments from the fits are shown in Figure [| as solid lines for fits to MLLA-0 spectra and as dashed 
lines for fits to the MLLA-M spectra. The data points are from our fit to e + e~ data [|[ |o). 

It may be concluded from Figure ^ that the MLLA-0 model gives a good description of the e + e~ data for 
the mean and skewness. There are however discrepancies in a and kurtosis as is also evident in Figure a is 
smaller, and the model is more platykurtic than the data. The MLLA-M model gives a poorer description of 
all variables at low Eqm but approaches both the MLLA-0 predictions and the data at high Ecm- 

This is contrary to the conclusions reported by Lupia and Ochs |fllf . In that paper the moments of the 
experimental distributions corrected for the mass effects were compared to the predictions for the limiting 
spectra and found to be in good agreement. It should be noted that in |li"| ] an extrapolation was made into 
regions where no experimental measurement exists to allow the moments to be calculated. In addition, the 
agreement of the MLLA limiting spectra with the mass corrected experimental £ distribution is poor. The 
spectra at low £ are well described but the region around the peak and at higher £ is generally not well 
reproduced. The studies reported here avoid the need to extrapolate into regions not experimentally measured. 

LPHD postulates that a constant factor relates the MLLA predicted spectrum to the experimental spectrum, 
independent of the Ecm- In the spectra discussed above, the maxima of the MLLA predictions have been scaled 
to the data. This scale factor, for the MLLA-0 model, is shown in Figure [| as a function of Ecm- It falls from 
w 1.4 at the lowest Ecm to ~ 1.2 at the highest. A similar result is obtained for an area normalisation over 
approximately three units of £ p for the MLLA-M case. Thus our observations disagree with the LPHD postulate 
of a constant scale factor. This suggests that at the currently accessible experimental energies the theoretical 
spectra are still influenced by the Qq cut-off. 

Similar comparisons of e + e~ data with MLLA predictions have been made by the DEPHI |l2| and OPAL |L0| 
collaborations. CDF has examined jets produced in pp interactions |T^ |. In contrast to the analysis presented 
here, no explicit analysis of the higher moments of the £ distribution has been made. However, where the 
analyses can be compared, the results are in agreement with those presented here. In particular, A is close to 
250MeV, the shape of the MLLA predicted spectrum is close to that of the data but disagrees in detail, and 
the value of the LPHD constant changes with energy. 

5 Conclusions 

The MLLA predictions of the logarithmic scaled energy spectra, and in particular their moments, have been 
investigated as function of energy scale. Various theoretical approaches have been compared and differences 
discussed. We note that care should be taken in the comparisons to ensure that a consistent approach is 
maintained; in particular attention should be paid to the range of application in £. There is a large discrepancy 
in the skewness and kurtosis at low Q but the various predictions converge as the asymptotic limit is approached. 



The theoretical results are compared to measurements taken in c + c~ annihilation experiments. The calcu- 
lations are generally in good agreement with the data. It is observed that the limiting spectra is preferred over 
the predictions of the truncated cascade (Qo ^ A.) The introduction of the mass term has a large effect at all 
but the highest Ecu ■ The data is consistently broader than the limiting spectra over the energy range studied 
here. 

The normalisation factor of LPHD between the data and the MLLA theoretical predictions is not constant. 
Contrary to expectation it displays a dependence on Eqm, decreasing as Eqm increases, thus suggesting a 
residual influence of the Qo cut-off. 
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Figure 1: The logarithmic scaled energy (a)-(c) and scaled momentum (d) spectra of partons as calculated in the 
MLLA formalism. Unless otherwise stated values of A = 250 MeV and Qo = A (the so-called limiting spectra) 
at Q — 91.2/2 GeV were used to calculate the curves, (a) The three curves are predictions at three different 
energy scales, Q = 183.0/2, 91.2/2 and 14.8/2 GeV represented by the full, dashed and dotted lines respectively, 
(b) The two curves are predictions for two different values of Qo — A and 2A represented by the full and dashed 
lines respectively, (c) The three curves are predictions for three different values of A = 50, 250 and 400 MeV 
represented by the full, dashed and clotted lines respectively. (d)The logarithmic scaled momentum spectra of 
partons as calculated in the MLLA formalism. The two curves are predictions for a massless (full line) and 
massive assumption (dashed line). 




Figure 2: The evolution of the mean, a, skewness and kurtosis as a function of Q for A = 250 McV. The dashed 
and dash-dotted lines are the analytic calculation of Fong and Webber for the gluon and quark respectively, the 
dotted line from the analytic calculation of Dokshitzcr et ai, and the full line from fitting a distorted Gaussian 
around ±lcr of the mean of the MLLA-0 spectra. 




Figure 3: e+e - data compared with MLLA-0 spectra ( dashed histogram), MLLA-M spectra( solid histogram 
) and MLLA Q Q = 2 A with A = 50 MeV ( thin solid line). The thick solid line is from the fit of a distorted 
Gaussian that is used to determine the moments for the data. 




Figure 4: Moments as a function of Ecm{= 2Q). The points with error bars are from fits of a distorted Gaussian 
to e+c~ data; circles: TASSO; square: TOPAZ; triangles: OPAL. The solid lines are moments from a fit of a 
distorted Gaussian to the MLLA-0 spectrum. The dashed lines are moments from a fit of a distorted Gaussian 
to a MLLA-M spectrum. 
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Figure 5: Scale factor as a function of Ecm(= 2Q). The full symbols are from normalisation to the same 
maximum height for data and MLLA-0 spectra. The open symbols (displaced by 1 GeV for clarity) are from a 
normalisation of the areas of data and MLLA-M spectra; circles: TASSO; square: TOPAZ; triangles: OPAL. 
The line is a linear fit of the height-based scale factor to log(-E'cM)- 



